Sarma Phase in Trapped Unbalanced Fermi Gases 
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We consider a trapped unbalanced Fermi gas at nonzero temperatures where the superfluid Sarma 
phase is stable. We determine in particular the phase boundaries between the superfluid, normal, 
and phase-separated regions of the trapped unbalanced Fermi mixture. We show that the physics 
of the Sarma phase is sufficient to understand the recent observations of Zwierlein et al. [Science 
311, 492 (2006); Nature 442, 54 (2006)] and indicate how the apparent contradictions between this 
experiment and the experiment of Partridge et al. [Science 311, 503 (2006)] may be resolved. 
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Introduction. — In the last two years impressive ex- 
perimental progress has been made in the field of ul- 
tracold Fermi gases, in part due to the possibility to 
tune the interatomic interaction strength by means of 
a Feshbach resonance. This has led to the study of the 
crossover between the Bose-Einstein condensation of di- 
atomic molecules and the Bose-Einstein condensation of 
atomic Cooper pairs, the so-called BEC-BCS crossover 
P, 13, Sj • Most recently, two experimental groups have 
gone even a step further by obtaining also full control 
over the polarization of the Fermi mixture. This has al- 
lowed for the study of degenerate Fermi gases with imbal- 
anced spin populations, which is a topic of great interest 
in many areas of physics ranging from condensed-matter 
physics to nuclear and astroparticle physics. These pi- 
oneering experiments by Zwierlein et al. 0, Q and Par- 
tridge et al. [p induced a flurry of theoretical activity 

0, i i III Q E E E IS, E H E I . 

The ^Li experiments of Zwierlein et al. and Partridge et 
al. both revealed that superfluidity in an ultracold Fermi 
gas is maintained upon going to an unequal mixture of 
two spin states. However, rather contradictory results 
are obtained by both experimental groups for the behav- 
ior of the Fermi mixture as a function of the population 
imbalance. Zwierlein et al. observe a phase transition 
between the superfluid phase and the normal phase at 
a high critical polarization of about 70%, whereas Par- 
tridge et al. seem to observe a transition between two 
different superfluid phases at a low critical polarization 
of about 10%. One important aim of the present Letter 
is, therefore, to propose a single theoretical picture in 
which the qualitative differences in the observations by 
the two experimental groups can be understood as two 
different sides of the same coin. Based on this picture, 
we also make a detailed quantitative comparison with the 
experiments of Zwierlein et al. 

The main results of our mean-field calculations are pre- 
sented in Fig. 1. Here we show the universal phase dia- 
gram of a trapped Fermi gas in the unitarity limit as a 
function of temperature and polarization. This phase di- 
agram is universal in the sense that it does not depend on 
the total number of fermions or the trap geometry. Note 
that in determining the phase diagram we have neglected 



Normal 
Phase 



0.75 



0.5 ■ 



0.25 




FIG. 1: (Color online) Universal phase diagram of a trapped 
unbalanced Fermi gas in the unitarity limit. The polarization 
P is given by (A'"+ — N-)/{N+ + N-), where N± designates 
the number of fermions in each hyperfine state of the Fermi 
mixture. The temperature T is scaled with the critical tem- 
perature Tc of the balanced Fermi gas. The solid horizontal 
line gives the temperature that is used for the comparison 
with the experiments of Zwierlein et al. presented in Figs. 2 
and 3. The points a, b, and c correspond to the polarizations 
used in the density profiles of Figs. 2a, 2b, and 2c. 

fiuctuations, which are known to be quantitatively im- 
portant in the unitarity limit. However, fluctuations are 
not expected to alter the topology of the phase diagram 
in this case. Fig. 1 reveals that there is a tricritical point 
in the trapped Fermi mixture, which is well-known for 
the homogeneous case [2^ 22^ i2M ' but has to the best of 
our knowledge not been studied yet for the harmonically 
trapped situation. In the normal phase the gas is in its 
normal state throughout the trap. In the Sarma phase 
the Fermi gas has a shell structure, in which the core of 
the trapped gas is superfluid, whereas the outer region 
is normal. Furthermore, the normal-to-superfluid transi- 
tion as a function of the position in the trap is of second 
order. Since the superfluid order parameter A vanishes 
continuously at the transition, we have for nonzero po- 
larizations always a region in the trap where |A| is so 
small that it results in a gapless superfluid with nega- 
tive quasiparticle excitation energies, as flrst studied by 
Sarma [24]. Since |A| increases towards the center of the 
trap, it is also possible for small polarizations and low 
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temperatures that the superfluid becomes gapped in the 
center of the trap. This leads to a gapped BCS super- 
fluid core with a gapless Sarma superfluid and normal 
shell surrounding it 2j]. Finally, in the phase-separated 
region of the phase diagram, the superfluid core and the 
normal shell of the gas are separated by a first-order tran- 
sition as a function of position, which implies that A goes 
discontinuously to zero at a certain cquipotcntial surface 
in the trap. 

Fig. 1 allows for a natural explanation of the qualita- 
tive differences in the observations by the two experimen- 
tal groups. More precisely, we will argue in the second 
part of this Letter that the experiments of Zwierlein et 
al. have observed the transition from the normal phase to 
the Sarma phase, implying that these experiments have 
been performed above the temperature of the tricriti- 
cal point. Moreover, we suggest that the experiments 
of Partridge et al. have been performed in the tempera- 
ture regime below the tricritical point, since these exper- 
iments appear to see the transition between a non-phase- 
separated and a phase-separated superfluid phase. Note 
that this explanation is fundamentally different from 
the proposal of Machida et al, who obtain a different 
phase diagram by considering the Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) phase, which has a spatially vary- 
ing superfluid order parameter, rather than phase separa- 
tion [16(] . As a result they obtain a Lifshitz point instead 
of a tricritical point. They argue that both experimen- 
tal groups operate beneath this Lifshitz point, making 
it difficult to explain the qualitative differences observed 
by both groups. Moreover, there is presently no exper- 
imental indication for the presence of the FFLO phase. 
Finally, we would like to stress, that the local-density ap- 
proximation that was used to obtain the phase diagram 
in Fig. 1 is not sufficiently accurate to describe all aspects 
of the experiments of Partridge et al. 13, l3, i3| ■ For this 
reason we consider here from now on only the experiment 
of Zwierlein et al. , where a quantitative comparison with 
our theory turns out to be possible. 

Universal phase diagram. — To obtain the phase dia- 
gram of Fig. 1 we use the mean-field theory for the Sarma 
phase in the local-density approximation as described by 
Houbiers et al. [25|. There the mean-field theory was 
applied to superfluid ^Li in the BCS-limit. Here we in- 
corporated the relevant physics of the unitarity limit by 
using a generalization of the approach put forward by 
Fregoso and Baym [11]. As a result, we end up with the 
following thermodynamic potential 



^ ( Ek - - TltJk 



k,o- 



2ek 



>-2 



(1) 



where the atomic dispersion is Ck = h^\i^ /2m with m the 
fermionic mass, A is the BCS order parameter, cr = ± 



denotes the two hyperfine states of the Fermi mixture, 
N„ is the number of atoms in each hyperflne state, 
and (3 = l/fc^T. Moreover, the average renormalized 
chemical potential /i' is given by fj! = -\- /^+)/2 
and the dispersions fiuj^^ u of the Bogoliubov quasipar- 
ticles satisfy ?iWk,o- = — A*'f)/2 -I- ?iwk with Tilo\^ = 

v/(ek - fJ-'y + |Ap. Finally, the densities n„ = N„/V 
are given by 

na^^Y. {l^k|'iV(?ia;k,.) + bkHl - A^(?iwk,^.)]} , 

k 
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where V is the volume of the mixture and N{hu}k^a-) = 
l/[exp(/3?ia;k,<T) + 1] are the Fermi distributions for the 
Bogoliubov quasiparticles. The BCS coherence factors 
|ukp and |wkP are determined by the relations |wkP = 
[1 -t- (ek - ^l')/huJk]/2 and \uk\^ + \vk\^ = 1. 

The renormalized chemical potentials are given by 
fi'^ = fj,^ — hT,^, where fj,^ are the chemical potentials 
and the selfenergies are approximated at unitarity by 



(p - PBcsjy — — "-<T • 



TO 



2fi' 



(3) 



This expression can be understood as follows. Since the 
average kinetic energy of the atoms involved in the inter- 
action is given by 2/x', the effective interaction strength 
between the atoms is expected to be proportional to 
l/\/2fi' [2^. To understand also the proportionality 
constant we make use of the fact that in the case of 
equal density and at zero temperature, Eq. (2) results 
in /i'^ = (1 -f f3Bcs)^F, where we introduced the Fermi 
energy e^? of a balanced Fermi gas and we also recall 
that in the BCS theory Pbcs - -0.48 Substituting 
this result for n'^ into Eq. (3), we obtain for the chemi- 
cal potential = l^'a + KTia = (1 + /?)eF, which agrees 
with the exact result from Monte-Carlo calculations for 
P ~ -0.59 [27, 



Taking the derivative of the thermodynamic potential 
with respect to A* and equating it to zero gives us the 
BCS gap equation 
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Note that in deriving this gap equation we did not dif- 
ferentiate the last term in the right-hand side of Eq. ([T]) . 
Differentiating also this term results in ffuctuation cor- 
rections to the mean-fleld theory on which we comment 
at the end of this Letter. 

The above theory is valid for a homogenous Fermi 
mixture in the unitarity limit. To account for the pres- 
ence of an axially symmetric trapping potential V{'x) = 
TO(wj_r2 -I- ujzz'^)/2 we use the local-density approxima- 
tion. This means that the theory is locally homogeneous 
with a spatially varying chemical potential. As a re- 
sult, the local renormalized chemical potential is given by 
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FIG. 2; Typical density profiles for trapped unbalanced Fermi gas above the tricritical point. The axial and radial directions 
of the trap are scaled in such a manner that the trap becomes effectively isotropic with a frequency lo — {Ld^uOzY^^ and a size 
I = {Ti/mujY^^ . The temperature and the polarizations used for the three panels are indicated in Fig. 1. The total number of 
atoms is A'^ = 1.4 ■ 10^. The dashed lines in panel (c) show the ideal gas results. 



/i^ (x) = — ^(x) — fiT,^ (x) . For a balanced Fermi gas at 
zero temperature, we thus have /io- — T^(x) = (l+/3)eF(x) 
and therefore we retrieve, as desired, the exact density 
profile in this case. Note that in the outer region of the 
gas cloud the renormalized chemical potentials become 
negative, so that 2/x'(x) in Eq. (3) is no longer a good 
measure for the avarage kinetic energy of the interacting 
fermions. Therefore, we then take SksT instead of 2fi' 
as an appropriate measure for the kinetic energy. 

We are now in the position to explain how we obtain 
the phase diagram of Fig. 1 . We first determine the line 
between the normal and the two superfluid phases. This 
is achieved by solving the BCS gap equation in the cen- 
ter of the trap and finding the temperature at which the 
BCS order parameter vanishes. Inspection of the thermo- 
dynamic potential reveals that the vanishing of the order 
parameter can occur continuously or discontinuously, i.e., 
by a second-order or a first-order phase transition. If the 
transition is of first/second order, we go from the normal 
to the phase-separated/Sarma phase. At the triciritical 
point these two different kinds of transitions merge. 

So far, we only looked at the center of the trap, but 
the tricritical condition can also be satisfied at a certain 
equipotential surface outside the center of the trap. This 
gives us a point on the Sarma-to-phase-separation line. 
To see this, consider a point on this line and raise the 
temperature slightly. This changes the tricritical tran- 
sition outside the center of the trap into a second-order 
transition slightly closer to the center of the trap, which 
means that the gas is in the Sarma phase. In a similar 
way, a slightly lower temperature leads to a first-order 
transition as a function of position in the trap, i.e., to 
the phase-separated phase. 

Comparison with experiment. — We now compare our 
theory with the experiments of Zwierlein et al. First, 
we show in Fig. 2 three typical density profiles of the 
gas: Two in the superfluid Sarma phase and one in the 
normal phase. Using different methods, similar density 



profiles have recently also been obtained by Yi and Duan 
im and Chien et al. The most striking feature 

in the density profiles is the 'bulge' in the minority and 
majority profiles in the Sarma phase. This 'bulge' is a 
direct consequence of the presence of the condensate of 
Cooper-pairs and was indeed one of the most important 
findings of the experiment. It shows that for an unbal- 
anced unitarity gas, in contrast to an unbalanced Fermi 
gas in the BCS limit [IBj, the condensate of Cooper-pairs 
has a very strong effect on the atomic density profiles. In 
the normal phase the attractive effects of the selfenergies 
given in Eq. (3) also slightly increase the central densi- 
ties of the two species. However, this occurs always in a 
smooth featureless manner whereas in the experiments of 




FIG. 3: (Color online) The radial size R of the Cooper-pair 
condensate (solid), and the minority (dashed) and majority 
(dot-dashed) gas clouds as a function of polarization for the 
temperature shown in Fig. 1. All radial sizes are scaled with 
the radial Thomas-Fermi radius of the balanced Fermi gas 
Rtf. Also shown is the experimental data of Zwierlein et al. 
.5]. For the radial size of the majority cloud the black points 
assume a hydrodynamic expansion, whereas the grey points 
assume a ballistic expansion. 
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Zwierlein et al. a distinct feature is seen in the majority 
density profile at the edge of the minority cloud. This is 
maybe due to the fact that our mean-field-like selfener- 
gies are intended to exactly incorporate the (3 parameter 
of the balanced unitarity gas, but only approximately ac- 
count for the strong correlations in the normal phase of 
the unbalanced unitarity gas. Another explanation could 
be that in the experiment the density profiles are deter- 
mined after an expansion. This affects in particular the 
majority density profile, since the outside of the majority 
cloud expands ballistically, whereas the inner part inside 
the minority cloud expands hydrodynamically and thus 
faster than the outer part. The expansion is therefore 
not determined by a single scale factor and this may lead 
to a pile up of atoms in the transition region. 

Besides the above qualitative comparison, we can also 
make a quantitative comparison with the experiment of 
Zwierlein et al. by determining the radial size of the 
Cooper-pair condensate and the radial sizes of the mi- 
nority and majority gas clouds. The radial size of the 
Cooper-pair condensate follows directly from the point 
where A(x) = 0, but the determination of the minority 
and majority radii is somewhat more complicated be- 
cause we are working at a nonzero temperature and the 
density profiles always have a gaussian tail. For sim- 
plicity we determine these radial sizes by the conditions 
/ijj(x) = 0, which give the correct results at T = 0, but 
underestimate the radial size of the density profiles at 
nonzero temperatures. The results of this procedure to- 
gether with the experimental data are shown in Fig. 3. 
In general, the agreement with experiment is very good, 
which confirms our picture that the experiment is oper- 
ating above the tricritical point. As expected, at small 
polarizations the radial size of the majority cloud agrees 
best with the experimental data obtained by assuming a 
hydrodynamic expansion, whereas for large polarizations 
our results approach the experimental data obtained by 
assuming a ballistic expansion. 

Discussion. — The main remaining problem of our 
theory is the high absolute value of the temperatures 
above the tricritical point. The calculated temperatures 
are typically a factor of 5 higher than what is found in 
the experiments. An important consequence of the high 
temperature is that at the Sarma-to-normal transition 
the densities of the two spin states are not equal in the 
center of the trap in contrast to what is found in exper- 
iments. The high absolute value of the temperature is 
the result of neglecting fluctuations, that substantially 
shift down the tricritical temperature. Including fluctu- 
ations in the case of a balanced Fermi gas reduces the 
critical temperature in the unitarity limit by a factor of 
3, and fluctuations are expected to be even more impor- 
tant in the unbalanced case. One reason for the latter is 
that the Sarma phase is a gapless superfluid, whereas the 
BCS phase at P = has a gap. However, theoretically 
the study of fluctuation effects is rather challenging for 



an unbalanced Fermi mixture at unitarity, since the usual 
Nozieres-Schmitt-Rin k ap proach Q has some unphysical 
features in this case [20|. Therefore we are developing 
a more advance theory including fluctuations. This will 
determine the location of the tricritical point more accu- 
rately, but, as mentioned previously, is not expected to 
alter the topology of the phase diagram. 
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